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PACS 7 1 . 3 . +h - Metal insulator transitions and other electronic transitions 

PACS 7 1 . 1 . HF - Non-Fermi-liquid ground states, electron phase diagrams and phase transitions in model 
systems 

PACS 7 1 . 3 8 . -k - Polarons and electron-phonon interactions 

Abstract. - We re-examine the Peierls insulator to Mott insulator transition scenario in the one-dimensional 
Holstein-Hubbard model where, at half-filling, electron-phonon and electron-electron interactions compete 
for establishing charge- and spin-density-wave states, respectively. By means of large-scale density-matrix 
renormalization group calculations we determine the spin, single-particle and two-particle excitation gaps 
and prove — in the course of a careful finite-size scaling analysis — recent claims for an intervening metallic 
phase in the weak-coupling regime. We show that for large phonon frequencies the metallic region is even 
more extended than previously expected, and subdivided into ordinary Luttinger liquid and bipolaronic liquid 
phases. 



The challenge of understanding the subtle interplay of 
electron-electron and electron-phonon interaction effects in 
low-dimensional condensed matter systems, such as conju- 
gated polymers, charge transfer salts, inorganic spin-Peierls 
compounds, halogen-bridged transition metal complexes, fer- 
roelectric perovskites, or organic superconductors, [1-4] has 
stimulated intense work on generic fermion/spin-boson models. 
In this respect the one-dimensional (ID) Holstein 1 -Hubbard 2 
model (HHM) is particularly rewarding to study. [5-14] It ac- 
counts for a tight-binding electron band, an intra-site Coulomb 
repulsion between electrons of opposite spin, a local coupling 
of the charge carriers to optical phonons, and the energy of the 
phonon subsystem in harmonic approximation: 



In Eq. (1), n i(J = c^ ia c ia , where c\ a (c ia ) creates (annihilates) a 
spin-CT electron at Wannier site i of an ID lattice with N sites, 
and b\ are the corresponding bosonic operators for a dis- 
persionless phonon with frequency qjq. 

The physics of the HHM is governed by three compet- 
ing effects: the itinerancy of the electrons (oc t), their on- 
site Coulomb repulsion (oc U), and the local electron-phonon 
(EP) coupling (oc e p ). Since the EP interaction is retarded, the 
phonon frequency (u>q) defines a further relevant energy scale. 
Hence one is advised to introduce besides the adiabaticity ratio, 



a = ujo/t , 
two dimensionless coupling constants: 
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'The Holstein model (cf. T. Holstein, Ann. Phys. (N.Y.) 8, 325 (1959); 
ibid. 8, 343 (1959)) has been studied extensively as a paradigmatic model for 
polaron formation in the low-density limit. For commensurate band fillings the 
coupling to the lattice supports charge ordering. 

2 The Hubbard model (cf. J. Hubbard, Proc. Roy. Soc. London, Ser. A 
276, 238 (1963)), originally designed to describe ferromagnetism of transition 
metals, has more recently been used as the probably most simple model to 
account for strong Coulomb correlation effects in a great variety of materials. 



u = U/4t, and g 2 = e p /uj or \ = e p /2t. (3) 

Both Holstein and Hubbard interactions tend to immobilize 
the charge carriers, and even may drive a metal-insulator tran- 
sition at commensurate band fillings. For the half-filled band 
case (one electron per lattice site), most previous analytical 
and numerical studies of the HHM reveal that Peierls insu- 
lator (PI) or Mott insulator (MI) states are favored over the 
metallic state at zero temperature. Whereas the PI is charac- 
terized by a distortion of the underlying lattice accompanied by 
dominant charge-density-wave (CDW) correlations, the Mott 
insulator is basically a spin-density-wave (SDW) state without 
any lattice dimerization. The physical excitations differ accord- 
ingly: while "normal" electron-hole excitations are expected in 
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Fig. 1 : Schematic representation of the Peierls insulator to Mott insu- 
lator quantum phase transition at fixed A and a in the strong coupling 
small-to-intermediate phonon-frequency regime of the ID half-filled 
Holstein-Hubbard model. The PI exhibits long-range order (LRO) be- 
cause a discrete symmetry is broken; the SDW MI has continuous 
symmetry, i.e. there is no LRO in ID. A s (dashed line) and A C1 (solid 
line) denote the spin and charge excitation gap, respectively. For finite 
periodic chains with N = An, the transition could be identified by 
a ground-state level crossing associated with a change in the parity 
eigenvalue P = rtl, where the site inversion symmetry operator P is 
defined by Pc\ a P^ = cj v _ l+1(T (cf. Refs. [5, 16]). 

the PI phase, charge (spin) excitations are known to be mas- 
sive (gapless) in the MI state, at least for the ID Hubbard-only 
model. Thus, varying the control parameter u/X, a cross-over 
from standard quasi-particle behavior to spin-charge separation 
might be observed in the more general HHM (see Fig. 1). This 
scenario has been corroborated by extensive Lanczos diago- 
nalization and density-matrix renormalization group (DMRG) 
studies. [5,7, 15] In particular the existence of a single quantum 
critical point, separating Peierls- and Mott-insulating phases 
at u/X ~ 1, has been confirmed by up-to-date stochastic 
series expansion (SSE) quantum Monte Carlo (QMC) calcu- 
lations. [12] But this only holds in the strong-coupling and 
adiabatic-to-intermediate phonon-frequency regime. 

For the pure (spinless and spinful) Holstein model it is 
well known that quantum phonon fluctuations may destroy 
the Peierls phase, provided the EP interaction is not too 
strong. [17-22] In this case, a Luttinger liquid phase exists be- 
low a critical coupling g c , where X c — > as a — > 0. 3 For 
this reason, it is natural to examine the stability of such a 
metallic state as the Coulomb interaction is turned on. (We 
call metal any phase with gapless charge excitations.) An 
intervening metallic phase would of course prevent a direct 
insulator-to-insulator transition. Recent numerical investiga- 
tions of the ID HHM, based on variable-displacement Lang- 
Firsov [6], DMRG [9, 11], and SSE QMC [8, 12] approaches, 
give strong evidence that the CDW-SDW transition does in- 
deed split into two subsequent CDW-metal and metal-SDW 
transitions in the weak-coupling intermediate-to-large phonon- 
frequency regime. 4 To map out the phase diagram is techni- 
cally challenging, however, since the energy scales are not well 

3 The metal-insulator transition at g c is expected to be of Kosterlitz- 
Thouless type, at least for u>a — > oo (cf. Refs. [23], and [24]). 

4 A intermediate metallic phase has been shown to exist also for the D = oo 
HHM (see, e.g., Ref. [25] and references therein). 
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Fig. 2: (Color online) Finite-size scaling of spin (left and right tri- 
angles; red dashed lines) and charge (diamonds, triangles up, circles; 
blue solid lines) excitation gaps in the ID HHM with OBC. Filled 
(open) symbols belong to insulating (metallic) states, see Fig. 4. 

separated in this region. Furthermore, it is not obvious which 
physical quantities are most suitable to distinguish the various 
phases. In this respect, local magnetic moment and effective 
electronic hopping integral [6], charge-, spin- and (supercon- 
ducting) pairing-correlation functions [9], or Luttinger liquid 
parameter [8] and charge/spin susceptibilities [12] have been 
proposed. 

In this work our analysis of the weak-to-intermediate cou- 
pling regime of the ID HHM is based on the behavior of var- 
ious spin and charge excitation gaps. These many-body gaps 
can be calculated by DMRG with high precision on a sequence 
of sufficiently large lattices. This allows for a reliable finite-size 
scaling of the gap data using their expected asymptotic behav- 
ior A(N) = A(oo) + b/N + c/N 2 for N > 1. A detailed 
description of our fermion-boson DMRG pseudo-site method 
(which treats electron and phonon degrees of freedom on an 
equal footing in the whole parameter regime) can be found in 
Ref. [26]. In the numerical work we use typically 3-5 bosonic 
pseudosites and employ open boundary conditions (OBC). 

The charge and spin excitation gaps are determined from 

A Cl = E+(l/2) + £ -(-l/2) - 2E (0) (4) 

and 

A s = £ o (l)-£o(0), (5) 

respectively, where Eq ± \s z ) is the ground-state energy at 
(away from) half-filling with N e — N (N e = N ± 1) particles 
in the sector with total spin-z component S z . In the various 
phases expected to occur in the ID HHM A s corresponds to 
the energy of the lowest spin excitation (spin gap). The one- 
particle gap A C1 yields the energy of the lowest charge excita- 
tion (charge gap) in the Peierls band insulator, Luttinger liquid, 
and Mott insulator phases and is an upper limit for this charge 
gap in any other phase. Note that since we compare ground- 
state energies when calculating the charge and spin gaps, lat- 
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Fig. 3: (Color online) Finite-size scaling of the two-particle excitation 
gap A c , in the HHM. For comparison we include A C1 in the bipolaron 
liquid phase. For further explanation see text. 



tice relaxation effects arising from different particle numbers 
are included. 

Searching for metallicity in the HHM, we first consider the 
anti-adiabatic regime (a = 5), because large phonon frequen- 
cies will clearly act against any static Peierls ordering. Figure 2 
shows the scaling of A Cl and A s with system size N as the EP 
coupling A is lowered at fixed Coulomb interaction u = 0.25. 
For large A there is a CDW made up of "singlet bipolarons", i.e. 
a bipolaronic superlattice. Since the bipolarons that emerge are 
rather small objects, the finite-size dependencies of A C1 and 
A s are weak. This is just as in the strong-coupling non-too- 
anti-adiabatic regime [27] (cf. also Fig. 11 in Ref. [15]). With 
decreasing EP coupling, i.e. increasing ratio u/X, both A Cl and 
A s become smaller, acquire a notable finite-size dependency, 
and finally both scale to zero, indicating Luttinger-liquid metal- 
lic behavior. Since A Cl and A s converge to the same value as 
N — > oo, spin and charge degrees of freedom "stick together". 
Above a critical ratio u/X, however, a SDW (Mott insulating) 
phase is realized. In the MI we found a finite charge excitation 
gap, which in the limit X/u <C 1 scales to the charge gap of 
the Hubbard model, whereas the extrapolated spin gap remains 
zero (cf. the data for A = 0.05). This is in agreement with 
spin-charge separation. 

In order to gain deeper insight into the nature of the interme- 
diate phase between PI and MI we determine, besides A Cl , the 
two-particle excitation gap 



E 2 +(0) + E 2 ~(0)~2E (0) 



(6) 



This gap corresponds to the charge gap in a bipolaronic in- 
sulator and is an upper limit for it in any other phase. Of 
course, one- and two-particle excitation gaps should simulta- 
neously open if we enter the PI and MI phases. If the PI phase 
is a bipolaronic insulator (superlattice) rather than a traditional 
Peierls band insulator, mobile bipolarons may occur first in the 
dissolving process of the PI, as the X/u ratio is lowered. Such 
a bipolaronic metal/liquid phase will then be characterized by 
A C2 = but finite A Cl (and A s ). Adding/removing a sin- 
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Fig. 4: (Color online) Phase diagram of the half-filled Holstein- 
Hubbard model in the weak-coupling anti-adiabatic (a = 5) regime. 
Here filled squares, open triangles down, and open triangles up, de- 
note the PI, bipolaronic metal, and Luttinger-liquid metal phases, re- 
spectively. The extension of the two latter phases is marked by the 
intensely (weakly) shaded regions. Dashed (solid) lines designate the 
Peierls-intermediate state (intermediate phase-Mott) phase bound- 
aries obtained by Ref. [8]. The dotted line indicates u — A. The 
inset gives the (N — > oo) extrapolated values of the one-particle, two 



particle, and spin excitation gaps at u 
large as A Cl in the MI phase. 



1. Note that A c , is twice as 



gle particle from the metallic bipolaron phase is energetically 
costly because the bipolarons are (tightly) bound. A bipolaron 
as a whole, however, can be added or removed without effort. 

Figure 3 illustrates that this scenario holds in the anti- 
adiabatic weak-coupling regime. As the EP coupling gets 
weaker, we enter a region where A C1 > but A C2 — > (see, 
e.g., the data for A = 0.625). A C2 stays zero as A C1 vanishes at 
still smaller A (u fixed), until we enter the MI state. 

Figure 4 combines the findings of our scaling analysis for 
A s , A Cl and A C2 with the results of previous studies of the 
half-filled HHM based on susceptibility and Luttinger liquid 
parameter data [8, 12]. First, we note that the Pi-metal phase 
boundary is shifted compared to the results of Ref. [8], while 
the metal-MI transition line is the very same. That is, taking 
A C2 = as a criterion for the instability of the PI phase, we 
find an even larger region for the PI-MI intervening state. Sec- 
ond, within the metallic state, our data suggests a cross-over 
between a bipolaronic liquid ( A C2 = 0; A Cl , A s > 0) and a 



Luttinger liquid (A C2 = A C1 
The mean kinetic energy, 



0). 



t r \ 



(7) 



shown in Fig. 5 for u = 1, exhibits a pronounced minimum 
in the intermediate metallic phase [6]. The local magnetic mo- 



There is an ongoing discussion, whether superconducting pairing correla- 
tions may become dominant in the HHM between the CDW and SDW phases 
(cf., e.g. Ref. [9]). Here, we will not address the highly controversial issue of 
existence of superconducting ground states (off-diagonal long-range order) in 
a stricly ID electron-phonon model. 
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Fig. 5: Kinetic energy Ehm and local magnetic moment Lo at it = 1, 
a = 5. 



ment, 



(8) 



can be taken as as a measure for the "localization" of the elec- 
tron spin. For the pure half-filled Hubbard model it varies be- 
tween 3/8 (band limit) and 3/4 (atomic limit). As can be seen 
from Fig. 5, Lq is suppressed in the PI CDW phase because of 
the alternating arrangement of (almost) empty and double oc- 
cupied sites. Entering the metallic phase Lq is enhanced and 
reaches its maximum in the MI phase. Both Ekm and Lq are 
however not very suitable for fixing the phase boundaries. 

The corresponding results for the adiabatic regime (a = 0.5) 
are given in Figs. 6 and 7. Again we have strong evidence for 
an intermediate metallic state. The region where bound mo- 
bile charge carriers (bipolarons) exist, however, now is a small 
strip between the PI and metal phases only, and expected to 
vanish if the adiabaticity ratio a goes to zero. 6 If all energy 
scales — set by A, u, a — become small, the finite-size scaling 
is extremely delicate and it is difficult to resolve the differ- 
ent metallic phases. The main differences in comparison to 
the anti-adiabatic regime pertain to (i) the nature of the Peierls 
phase, which is a traditional band insulator in the adiabatic 
case, and (ii) the existence of a tri-critical point at moderate 
coupling strengths. 

To summarize, we complemented previous numerical inves- 
tigations of the half-filled Holstein-Hubbard model [6-9, 11,12] 
by a large-scale DMRG analysis of the weak electron-phonon 
and electron-electron interaction regime. By calculating spin 
and charge excitation gaps we confirmed the existence of an 
intermediate metallic phase in the cross-over region of the 
Peierls-insulator Mott-insulator transition. This phase is shown 
to be most extended in the anti-adiabatic limit of large phonon 
frequencies. If the PI typifies a bipolaronic superlattice, it gives 
way to a bipolaronic liquid composed of bound mobile (singlet) 
charge carriers as the EP coupling weakens. Reducing the EP 
to Coulomb interaction ratio further, a cross-over to a metal- 
lic state with unbound (weakly phonon-dressed) electrons (Lut- 
tinger liquid) takes place, until the MI state is finally reached 
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Fig. 6: (Color online) Finite-size scaling of the one- and two-particle 
excitation gaps in the adiabatic (a = 0.5) regime of the HHM. The 
inset gives the results in the PI phase A = 0.375. 




For the so-called adiabatic HHM with frozen phonons an additional bond- 
order-wave might occur within the intermediate phase. [28] 



Fig. 7: (Color online) Phase diagram of the HHM in the adiabatic 
weak-coupling regime. Again filled squares, open triangles down, 
and open triangles up, denote the PI, bipolaronic metal, and Luttinger- 
liquid metal phases, respectively. Dashed (solid) lines are taken from 
Ref. [8]; the dotted line gives u = A. Since the finite-size scal- 
ing becomes exceedingly costly in the vicinity of the tricritical point 
it ~ 0.25, A = 0.3375 given by Clay et al. [8], we cannot reliably 
resolve the phase structure there. 



at A < u. If the PI typifies a rather standard band insulator 
in the adiabatic, small phonon-frequency limit, the bipolaron 
liquid phase narrows substantially, but still a metallic phase ex- 
ists in between CDW and SDW states. Putting these findings 
together with the results previously obtained for the strong- 
coupling case, where a direct first-order PI-MI transition takes 
place, a rather complete picture of the physical properties of the 
ID HHM emerges. 
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